The three-body Coulomb scattering problem in discrete Hilbert-space basis 
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For solving the 2 — > 2, 3 three-body Coulomb scattering problem the Faddeev-Merkuriev integral 
equations in discrete Hilbert-space basis representation are considered. It is shown that as far as 
scattering amplitudes are considered the error caused by truncating the basis can be made arbi- 
trarily small. By this truncation also the Coulomb Green's operator is confined onto the two-body 
sector of the three-body configuration space and in leading order can be constructed with the help 
of convolution integrals of two-body Green's operators. For performing the convolution integral 
an integration contour is proposed that is valid for all energies, including bound-state as well as 
scattering energies below and above the three-body breakup threshold. 
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I. INTRODUCTION 

The three-body Coulomb scattering problem is one 
of the most challenging long-standing problems of non- 
relativistic quantum mechanics. The source of the dif- 
ficulties is related to the long-range character of the 
Coulomb potential. In standard scattering theory it is 
supposed that asymptotically the particles move freely 
what is not the case if Coulombic interactions are 
involved. The problem has been formulated in the 
Faddeev-Merkuriev theory jlj] on a mathematically sound 
and elegant way. This theory has been developed trough 
integral equations with connected (compact) kernels and 
equivalently transformed into configuration-space differ- 
ential equations with asymptotic boundary conditions. 
However, in practical calculations, so far only the latter 
version of the theory has been considered, since the ex- 
tremely complicated structure of the Green's operators 
in the kernels of Faddeev-Merkuriev integral equations 
has not yet allowed a direct solution. 

Recently a novel method was proposed for treating 
the three-body Coulomb problem by solving integral 
equations in a discrete Hilbert-space basis representa- 
tion. In particular the Faddeev-Noble equations j^] 
were solved in the Coulomb-Sturmian basis represen- 
tation. The Coulomb Green's operators occurring in 
the Faddeev-Noble integral equations were connected, 
via Lippmann-Schwinger-type equations, to the channel- 
distorted Coulomb Green's operators, which can then 
be calculated by a convolution integral of two-body 
Coulomb Green's operators. The method was elabo- 
rated first for bound-state problems ||] with (repulsive) 
Coulomb plus nuclear potentials, then it was extended for 



calculating p — d scattering at energies below the breakup 
threshold Q. In these calculations excellent agreements 
with the results of other established approaches were 
found and also the efficiency and accuracy of this new 
method were demonstrated. In a more recent work 
atomic bound-state problems with attractive Coulomb 
interactions were considered among others it showed 
the surprising result that in this scheme a much smaller 
number of angular momentum channels were needed to 
reach convergence than in previous methods. However, 
the validity of the approach along the Faddeev-Noble is 
restricted. Without limitations of the energy range it is 
strictly valid only for repulsive Coulomb interactions. 

The Faddeev-Merkuriev theory is more general as it 
comprises all kind of dynamical situations and applies to 
repulsive as well as attractive Coulomb interactions. In 
this article we consider scattering processes which start 
from two-body asymptotic channels, i.e. 2 — > 2 and 
2^3 scattering. We adapt the discrete Hilbert-space 
basis representation to the Faddeev-Merkuriev integral 
equations for those processes. We provide justification 
for the validity of the expansion method for all energies 
below and above the breakup threshold and for repulsive 
as well as attractive long-range forces. 

In the following section we recapitulate the ba- 
sis elements of the integral-equation formulation of 
the Faddeev-Merkuriev scattering theory for long-range 
forces. In Sec. Ill we present our method for solving 
these equations and prove the validity of the involved 
expansions. We summarize the results and give some 
conclusion in Sec. IV. 
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II. FADDEEV-MERKURIEV INTEGRAL 
EQUATIONS 

We consider three particles with masses m a , mp, m 7 , 
where (a, /?, 7) are permutations of (1,2,3). By x a ,y a 
we denote the usual Jacobi coordinates, and we will also 
use the notation X = {x a ,y a } £ R 6 . The three-body 
kinetic-energy operator is given by the following sum of 
two-body kinetic energy operators 



We assume that the particles interact via (long-range) 
Coulomb-like potential 



v a (x a ) 



(1) 



where e denotes the electric charge. The short-range 
potential v s J(x a ) is supposed to be square integrable, 
v a{ x a) € L2CR 3 ), and consequently it should possess 
the asymptotic behavior 

V S a r (Xa)~0(\x a \- 3 / 2 - e ), £>0, 

as \x a \ —>■ 00. 

The Faddeev integral equations are the fundamental 
equations for three-body problems. They were originally 
formulated for short-range interactions g). In this case, 
after one iteration they possess connected kernels, and 
consequently they are effectively Fredholm-type integral 
equations of the second kind. The fact that the homoge- 
neous equations allow nontrivial solutions only for a dis- 
crete set of three-body bound-state energies guarantees 
the existence of unique solutions of the inhomogeneous 
equations at scattering energies. For long-range interac- 
tions, however, the integral kernels of the Faddeev equa- 
tions become singular and the nice properties above do 
no longer hold. Therefore one cannot simply plug in a 
Coulomb-like potential into the Faddeev equations. 

A natural way of remedy is to split somehow the 
Coulomb-like interaction into short-range and long-range 
terms 



.,(0 



and apply the Faddeev procedure only to the short-range 
parts. Here, the superscripts s and I indicate the short- 
and long-range attributes, respectively. The three-body 
Hamiltonian is defined by 



H = H° 



.(«) 



,(0 



XI) 



Merkuriev introduced the splitting in the three-body 
configuration space R 6 via a splitting function £ a , 

v^\X) = v a (x a )( a (x ai y a ), 



v£'(X) = v a (x a )[l - Ca(x ai y a )}. 



The function £ Q is defined such that it separates the 
asymptotic two-body sector f2 a from the rest of the three- 
body configuration space. On the region of Q a the split- 
ting function £ a asymptotically tends to 1 and on the 
complementary asymptotic region of the configuration 
space it tends to 0. Rigorously, f2 Q is defined as a part of 
the three-body configuration space where the condition 

\x a \ < o(l + \y a \/a) v ,with a > 0, < v < 1/2, (2) 

(s) 

is satisfied. So, in fl a the short-range part Va coin- 
cides with the original Coulomb-like potential v a and in 
the complementary region vanishes, whereas the opposite 
holds true for v$- 

Wc define the long-ranged and the channel long-ranged 
Hamiltonians as 



and 



respectively, together with their resolvent operators 



and 



GV>(z) = {z-H<P)- 



The Hamiltonian is constructed to be a positive op- 
erator with continuous spectrum a(H^) = [0, 00) 
This continuous spectrum is structureless, its properties 
are similar to the spectrum of the three-body kinetic en- 
ergy operator. As far as the resolvent is concerned it 
means that (z) is an analytic operator- valued func- 
tion in the complex z plane except for the cut coinciding 
with a(HV'). This is a very important property, it en- 
sures that the Faddeev procedure becomes an asymptotic 
filtering ^ and consequently guarantees the asymptotic 
orthogonality of the Faddeev components. In Noble's 
approach Q| the splitting is performed in the two-body 
configuration space, like in Eq. (Q), thus the above condi- 
tion, without any restriction to the energy range, can be 
fulfilled only in case of repulsive Coulomb interactions. 
Now, if we perform the Faddeev decomposition of the 
three-body wave function | ty) , 

|tf a > = G« w W|*>, 

it will play the role of an asymptotic filtering. It leads to 
the decomposition of the wave function into three com- 
ponents 

I*) = \i>a) + \il>fi) + |^y>, 

with each of them describing different asymptotic chan- 
nels. The resulting equations for the components are the 
set of Faddeev-Merkuriev integral equations M 
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i¥=<* 



(3) 



where is an eigenstate of Hi 1 ' with energy E 



f /3i 



Since we consider here only 2 — > 2, 3 processes de- 
scribes asymptotically bound-state pairs in the two-body 
subsystem (3 with binding energy Epi. Merkuriev has 
shown in Ref. |7|J|] that after a certain number of itera- 
tions Eqs. (^) are reduced to Fredholm integral equations 
of the second kind with compact kernels for all energies, 
including energies below (E < 0) and above (E > 0) the 
three-body breakup threshold. Thus all the nice prop- 
erties of the original Faddeev equations established for 
short-range interactions remain valid also for the case of 
Coulomb-like potentials. 

It is important to note here that the resolvent Ga of 
the channel long-ranged Hamiltonian Ha is an opera- 
tor which, in contrast to the short-range case, cannot 
be constructed only from two-body quantities. To deter- 
mine Ga uniquely one should start again from Faddeev- 
type resolvent equations, or from the triad of Lippmann- 
Schwinger equations |10|, which is, in fact, the adjoint 
representation of the Faddeev operator |ll|] . The Hamil- 
tonian Ha , however, has the advantageous property that 
it supports bound states only in the two-body subsys- 
tem a, and thus only one kind of asymptotic channel 
is possible, the a channel. For such systems one sin- 
gle Lippmann-Schwinger equation treating properly the 
long-range part of H a is sufficient for a unique solution 
Jfjj. The appropriate Lippmann-Schwinger equation was 
proposed by Merkuriev and has the form 



G^(z)=G a a s (z)+G a a s (z)V: s G<i\z), 



(4) 



with a three-body potential V as , which, as \X\ — > oo, 
decays faster than the Coulomb potential in all direc- 
tion of the three-body configuration space R 6 : V as ~ 

oa-xr 1 -'), e>o. 

The operators G° s and V£ s are complicated three- 
body operators. However, in the region Q a they take 
relatively simple forms. The potential V as is defined by 



V 



u a 



(0 



.CO 



(5) 



where the auxiliary two-body potential uf (y a ) must be 
taken in such a way that it has the asymptotic form 
e Q (e ( g + e 7 )/|y Q |, as \y a \ — > oo. In this case 
U a , as \y a \ — > oo, decays sufficiently fast, namely 



Ua\y a ) 



U c 



0(\y a 



-2+v\ 



(6) 



which, considering (|^), means square integrability. So, in 
the two-body region f2 Q the auxiliary long-range poten- 
tial u2' compensates asymptotically the long-range tail 
of 

In this region also the Green's operator G a s (z) takes 
a simple form, it coincides with the channel distorted 



Coulomb Green's operator G a (z). The operator G a (z) 
appears as a resolvent of the sum of two commuting two- 
body Hamiltonians 



hy a )-\ 



.(0 



therefore it 



G a (z) = (z-h 
where h Xa — h® + v a and h y 

can be constructed as a convolution integral of two-body 
Green's operators 

G%(z) = G a (z) = ± I dz'g Xa (z-z')g y Jz'), (7) 
2vri J c 

where g Xa (z) = [z — /i^) -1 and g Va (z) = (z — h Va ) . 
The contour C is to be taken counterclockwise around 
the continuous spectrum of h Va so that g Xa is analytic 
on the domain encircled by C. 

The spectra of Ga and of G^ s should be of similar na- 
ture. Therefore, if we want that G a replaces G° s on the 
domain of Q a it must possess an analogous spectrum. 
Since Gat has been constructed such that the branch 
points of its spectrum are related to the two-body bound 
states generated by v a (x a ) the auxiliary two-body poten- 
tial Ua\y a ) must be taken in such a way that it does not 
support any bound state. If e a (e@ + e 7 ) > 0, it is easy to 
fulfill this requirement, if e Q (e^ + e 7 ) < the infinitely 
many bound states should be projected out leading to a 
nonlocal uf . 

All these properties of Eq. (||) lead to the following 
asymptotic behavior of the integral 



Q (X,E)(^-\f) +Y,M*~)Qcd(v«,E)($®-\f), (8) 



where / e L2 . Here . 



is the j-th bound-state wave 



The state |$^) is 



function of the pair a with energy E, 

a scattering eigenstate of the Hamiltonian H$ describing 
in the final state a Coulomb distorted propagation of the 
three particles. The state describes the scattering 

of the third particle with respect to the bound pair. Its 
asymptotic form is given as a product of two-body wave 
functions 



where \Xa) is a scattering eigenstate of h Va . The 
Coulomb spherical wave functions possess the asymp- 
totic behavior Q (X,E) — \X\~~ 5 / 2 , as \X\ — > 00, and 
Qaj(y a ,E) ~ \ya\~ 1 , as \y a \ -> 00, respectively 0. From 
its construction it follows that the state \F^ S ) — G a s \f) 
possesses an asymptotic behavior similar to \F a ) of Eq. 

The asymptotics (H) leads immediately to an analo- 
gous asymptotic form for the solutions to the Faddeev- 
Mcrkuriev equations 
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1p a (X) ~ Spafoiixpfop* {y/3,E - Eft) + 

^ fiaj {Xg)Qaj (Va, E)A a j t pj + 
3 

Qo(X, E)A a o t pi , 
where the scattering amplitudes are defined by 

7^q 



and 



A a0 ,p i = ^($ i ll-\vi s ^). 



(9) 
(10) 

(11) 



Finally we should note that the term wi s ^|^ 7 ), 7 ^ a, 
appearing at the right hand side of Eq. (g) is a square 
integrable function in the three-body Hilbert space: 
Va^lipj) £ ^(R 6 ). In fact, according to (ph, we need to 



consider only terms like Va (f>-yj(x-y) and Va 'Qo{X, E). 
The function Va 4>"ij{ x i) IS square integrable due to the 

(s) 

fast decrease of the potential Va in the coordinate x a 
and of the two-body bound-state wave function </> 7 j(x 7 ) 

in the coordinate x-y. In the term v& Qo(X, E) the po- 

(s) 

tential Va contains the splitting function £ a , which re- 
stricts the domain of integration to tt a . On this domain, 
however, Qq(X, E) itself is square integrable if v < 1/2, 
therefore Va^Qo(X, E) is also square integrable in R 6 . 



III. SOLUTION OF THE INTEGRAL 
EQUATIONS 



Note that, in the three-body Hilbert space, three equiv- 
alent bases belonging to fragmentations a, (3, and 7 are 
possible. 

We start from the observation of the previous section 
that each term Va \ipf), a 7^ 7, in Eqs. (|J) is a square 
integrable function and consequently can be expanded in 
L-2(R 6 ) in the convergent series 



71,771=0 



I ram) a a (nm\v^\tp-y) 



JVi\ 



(13) 



In the last line of this equation we have truncated the 
series at some finite Ni and represented the remainder as 
{r^ 1 ) = (1 — l^ 1 )?;^ |^i 7 ). Due to the square integrability 
of Va ^Tp^} the latter has the zero limit 



lim ||;- 



Will - 



(14) 



where ||/|| = (/|/) 1/2 stands for the Hilbert -space norm. 
Now, the Faddeev-Merkuriev integral equations ([|) take 
the form 



+ G«\E + iO)\r^). 



(15) 



To proceed further let us consider 
the terms n (nm\vg^ appearing in the second term 
on r.h.s. of (|l5|). This matrix element can explicitly be 
expressed by the integral 



In this section we consider a discrete Hilbert space 
basis representation of the Faddeev-Merkuriev integral 
equations. We will show that the truncation of the basis 
at finite N leads to a negligible error in the scattering 
amplitudes if N tends to infinity. 



A. Discrete Hilbert space basis representation of the 
Faddeev-Merkuriev integral equations 

We suppose that the states {|ra)} together with their 
biorthogonal partner {\n)} form a basis in the Hilbert 
space L2(R 3 ) of the two-body relative motion. Since 
the three-body Hilbert space is a tensor product of two 
such two-body spaces, an appropriate basis is the direct 
product 

\nm) a = \n) a <8 |m) a , (n, m = 0, 1, 2, . . .). (12) 

The completeness of the basis in the three-body Hilbert 
space L 2 (R 6 ) can be expressed in the form 



1 = lim 

N- 



N 

, E 

n,m— 



\nm) a a nm 



lim 1 

N^oc 



N 



Q (nm|^ s) |-0 7 ) = / dx a dy a a {nm\X)v<£> (X)ifj^(X). 



(16) 

The product a (nm\X)va\x) in the integral is square 
integrable and again can be expanded in Z/2(R 6 ) using 
the basis (O) in fragmentation 7 



a (nm\X)v^(X) = 



(nm|i^|n'm') 7 7 (n'm'|X) 



= a (nm\vi s h^\X) + 
a (nm\vi s \l-&\X). 



(17) 



Again the expansion is terminated at some finite N2, but 
this is chosen such that for all n,m < Ni the remainder 
becomes arbitrarily small: 



,(nm|«M(l- <e; e > 0. 



(18) 



Substituting Eq. ( |17| ) into Eq. ( |15[ ) we arrive at the fol- 
lowing representation 
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\1> a ) = <WI< )+ > + G^(E + *0)1>« V 1?^ 7 > 



7 1^7/ 

7^q 



+ Gg>(£ + iO)[\r^) + \r^ N *)], (19) in E< 1- ® with amplitudes 



three-body wave functions l^j) and I^q), respectively. 
This involves a similar asymptotic form of as for if> a 



where |r^> = l^v^ £ (1 - l? 2 )!^)- 

7#" 

Below we will show that the last term in Eqs. ([19]) gives 
a negligible contribution to the scattering amplitudes if 
Ni and N 2 tend to infinity. Indeed, from Eqs. ([ll]) and 
( |l8| ) it follows that for any e > there exist iVj, and N 2 
such that for all N > max{iVi, N 2 } the inequality 

ll^+r^||< £ 

holds. The last term on the right-hand side of Eqs. jis| ) 
can now be considered as an additional driving term. 
Then the solution of Eqs. (|l9| ) can be represented as a 
sum of two terms 



|V«) = IV£> + K>, 



(20) 



which satisfy the Faddeev-Merkuriev integral equations 
with degenerate kernels of rank N 

K) = M*2 + > + Gi «( E + «£° £ ^K) ( 21 ) 

7^q 



\rZ) = g«\e- 



G«\E + iQ)l»vUT.l?\R N 



7 1 7 
7^a 



(22) 



<• = £«(!> 



)) 



and 



i& = £<*«ol(l 



+ 



)) 



From the unitarity relation of the three-body wave func- 
tions we get the inequality 



aj 



Ei^oi 2 



||r 



JV112 



< 3e 2 



This means that with a proper choice of N the moduli 
of the amplitudes can be made arbitrarily small and the 
neglect of \R&) in ( ^0|) causes only a negligible error in 
the scattering amplitudes. 

So far we have reduced the problem of solving Eqs. 
(||) to the solution of Eqs. (^l]) with a degenerate kernel 
for \4>a)- Its asymptotic behavior is exactly the same as 
of the original equations (|), only the amplitudes A a j,pi 
and A a j pi have to be replaced by 



A N 



=E<e 

7#" 



Due to the degeneracy of the kernels these equations are 
automatically compact. Actually, only the first one is 
needed for practical calculations because the solution of 
the second one can be made arbitrarily small by a proper 
choice of the parameter N. To show this explicitly we 
rewrite Eqs. (pin in differential form 



(£-^)K>-iM s) ViW 



/ j 7 la, 
7^a 



and express its solution via the resolvent of the Faddeev 
matrix operator fill] 



7 

The resolvent components G*L are defined as the solution 
of the Faddeev-Merkuriev integral equations 



which, as before, are automatically compact due to the 
degeneracy of the kernel. The same arguments that al- 
lowed to arrive at the asymptotics (0) lead to a similar 
asymptotic form for the kernel Gg a (X, X',E + i0). One 



0a 

CO- 



needs only to replace the states |3>^ ) and |$^q ) by the 



and 



A N 



£<*!' 

7^q 



(O-i 



respectively. So, in terms of amplitudes A and A N we 
can state that for any e > there exists Nq such that for 
every N > Nq the relation 



aj 



\A a j,i3i — A^jjjii 2 + yj \ A a0i 



Pi 



aN |2 , o 2 



holds, i.e. /3i and ^4^ ^ as N goes to infinity tend 
to A a j.j3i and Aqq,/^ respectively. 

In calculating A^ ^ and A^ ^ we can again perform 

an approximation. Recalling the fact that Va^\ipy) is 
square integrable, similarly to the approximations made 
before in formulas ( |l3| ) and (17), we perform approxima- 
tions to the scattering amplitudes 



aNM 



and 



aNM 



7#" 



(23) 



(24) 
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which, as M goes to infinity tend to ^ and A^ ^, 
respectively. Summarizing our results we can conclude 
that 



aNM |2 







as A" and M go to infinity. Thus the reduction of the 
Faddeev-Merkuriev integral equations to matrix equa- 
tions and the use of the resulting finite basis representa- 
tion of the Faddeev components in calculating scattering 
amplitudes via ( |25| ) and is justified, the error made 
by using finite N and M values in the truncation can be 
made arbitrarily small. 

In practical calculations all what we need is to solve 
the matrix equations for tp = a (nfn\'i/j a ) 



where G^} — a (nm\G^a\n'm') a , v_a] = Q (nm|u£ s ' ) |7T.'m') 7 

and $^ = a (jvm\^a) , and calculate the scattering am- 
plitudes by 



and 



A a o. 



0i 



/ j —aO —a-y^ 



B. Evaluation of matrix elements and 

— a — a 

(s) 

While the matrix elements of the potential Va can al- 
ways be calculated, at least numerically, we need further 
approximations to calculate the L2 matrix elements G^} 
and fW. 

For calculating G^ we make use of the Lippmann- 
Schwinger equation (Q), and we show that this integral 
equation for the matrix elements can be reduced to an 
integral equation with a degenerate kernel. Eq. (|J) for 
GjP takes the form 

G&O = (^\G a a s \ri7n') + (^\G a a s V° s G^\nhn>). (25) 

Let us consider a typical integral appearing in the sec- 
ond term of the previous equation 

(f\GZV?GW\f) = J dy^dxlF™{X")V™{X")FV{X"), 

(26) 

where /,/' e L 2 , and and F£ s are defined in (ph. 
We split the integration domain into two parts as 



dX" 



'r 6 



dX" 



dx" 



The leading term of the second integral, due to (N), for 
large £ values takes the form 



no 



dy y 
js. 



dx x 2 



cxp{i2y/E^x 2 + y 2 } 
(x 2 + y 2 )( 6+e )/ 2 ' 



and the asymptotic analysis, which involves integration 
by parts, results 



In other words, the integral ( |26D can be confined to 
the asymptotic two-body sector S1 Q of the configuration 
space. Thus we can use the explicit forms of G° s and V® s 
give n by (^) and (||) valid on this domain and instead of 
( P5| ) we can use 

= (nrnlG^n/m') + (^\G a U a G^\n i 7n'), 

where the restriction of integration on domain $7 Q is as- 
sumed implicitly. Now, due to the square-integrability 
of U a with respect to the variable y a (see Eq. (|^)) and 
because of Eq. (J|), we can write 

U a G<£\nWn>) = l%U a G£\nWn') + {l-l% )U a G^\vhn'). 
Repeating the arguments that lead us to Eqs. (|l^) we get 

G w = G + G U a G^ l) + 

(nrn|G Q (l - l^)U a G^\vhn') + 
(7M\G a l%U a {l-l*)G£\n^i'). 

By the same argumentation as before one can easily show 
that the contribution of the last two terms to the solution 
G^} can be made arbitrarily small by a proper choice 
of the expansion parameter N. Thus, we have reduced 
the problem of solving Eq. (|2^) to the purely algebraic 
problem 

G (z) = G + G U a G (l) . 

±La —a ' l=La— —a 

A similar equation holds true for the wave function ma- 
trix elements 

= $ 4- G U a <&® 

-£-a -£-a ~ ±La¥- —a ! 

where <& a = (n\4> a ) <8> (fn\xa^} are direct products of two- 
body wave function matrix elements. 

For evaluating the matrix elements G Q we can use the 
contour integral representation (0) 



GJz) = 



1 

2tti 



dz' a (n\g Xa (z-z')\n') a a (m\g ya (z')\m') a . 



The contour G is to be chosen so that it encircles the 
spectrum of g Va and avoids the singularities of g Xa . For 
bound-state energies the spectra of the two Green's oper- 
ators are well separated and this condition can be fulfilled 
easily [||. Also, for scattering energies below breakup 



G 



threshold, where the bound-state pole of g Xa can meet 
the continuous spectrum of g Va , the contour integration 
can still be performed [0] . For other positive-energy scat- 
tering problems, however, the continua overlap and there- 
fore these contours are not viable. 

Fortunately, there exists a contour that is valid for all 
z of physical interest. Besides positive real values of z 
that arise for scattering above the breakup threshold, 
this contour can deal with complex values of z having 
negative imaginary parts that are needed for resonant- 
state calculations. In this approach one must first shift 
the spectrum of g Xa by taking z = E + ie with positive e. 
By doing so, the two spectra become well separated and 
the spectrum of g Va can be encircled. Next the contour 
C is deformed analytically in such a way that the upper 
part descends to the second Riemann sheet of g Va , while 
the lower part of C can be turned away from the cut [see 
Fig. 0]. The contour still encircles the branch cut singu- 
larity of g Va , but in the limit e — > it now avoids the 
singularities of g Xa . Moreover, by continuing to negative 
values of e, the branch cut and pole singularities of g Xa 
move to the second Riemann sheet of g Va and, at the 
same time, the branch cut of g Va moves to the second 
Riemann sheet of g Xa . Thus, the mathematical condi- 
tions for the contour integral representation of G a (z) in 
Eq. (0) can be fulfilled for all energies. In this respect 
there is only a gradual difference between the calcula- 
tions for bound and resonant states as well as scattering 
below and above the breakup threshold. 

IV. SUMMARY 

In this paper we have shown that in the kernels of the 
Faddeev-Merkuriev integral equations for 2 — > 2,3 scat- 
tering the short-range potentials can be approximated 
with the use of a discrete Hilbert-space basis. Then, to 
solve for scattering amplitudes one needs only the ma- 
trix elements of the channel long-ranged Green's oper- 
ator between basis states. For evaluating these matrix 
elements we use specific Lippmann-Schwinger equation, 
and we show that in this particular case the potential 
term in the kernel can again be approximated on the 
same finite basis. This truncation involves a reduction 
of the asymptotic channel long-ranged Green's operator 
onto the two-body sector of the three-body configuration 
space, where it can be represented in leading order by a 
convolution integral of two-body Green's operators. For 
performing the convolution integral we propose a contour 
that applies for all energies of physical interest. 

In principle any discrete Hilbert-space basis would be 
admissible but we expect that a Coulomb-Sturmian (CS) 
basis |l3|] would be the most advantageous in practical 
calculations. The reason is that such a basis allows for 
an exact analytic calculation of the two-body Green's 
operators and of other two-body quantities [Of. As a 
result the contour integral can be performed also in the 



practice. The practicability of this approach has been 
demonstrated already in |5m). The present paper pro- 
vides the formalism for extending the previous approach 
to scattering problems above breakup energies and to re- 
pulsive as well as attractive Coulomb-like interactions. 
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FIG. 1. Analytic structure of g Xa (z — z') g ya (z) as a func- 
tion of z' with z = E + is, E > 0, e > 0. The contour C 
encircles the continuous spectrum of g Va and avoids the sin- 
gularities of g Xa . A part of it, which belongs to the second 
Riemann sheet of g Va , is drawn as a broken line. 
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